A Tight-Binding Approach to Overdamped Brownian Motion 
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We present a theoretical treatment of overdamped Brownian motion on a multidimensional tilted 
periodic potential that is analogous to the tight-binding model of quantum mechanics. In our 
approach we expand the continuous Smoluchowski equation in the localized Wannier states of the 
periodic potential to derive a discrete master equation. This master equation can be interpreted in 
terms of hopping within and between Bloch bands and for weak tilting and long times we show that 
a single-band description is valid. In the limit of deep potential wells, we derive a simple functional 
dependence of the hopping rates and the lowest band eigenvalues on the tilt. We also provide general 
expressions for the drift and diffusion in terms of the lowest band eigenvalues. 
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I. INTRODUCTION 

Brownian motion on a periodic potential has been used 
to describe a wide range of nonequilibrium transport phe- 
nomena where there are significant thermal fiuctuations 
[l,,^]. For example, this arises in Josephson junction ar- 
rays 3 , surface diffusion 0| , particle sorting experiments 
[B43]j ttis mechanochemical coupling of molecular 
motors (8l-[l3|. The standard theoretical description of 
these systems is based on a continuous diffusion equa- 
tion that governs the dynamics of the probability den- 
sity and can be used to determine the average of macro- 
scopic observables T] . In the limit of deep potential wells, 
it is physically intuitive that this continuous diffusion 
equation can be approximated by a discrete equation de- 
scribing infrequent hopping transitions between poten- 
tial wells. Discrete treatments have been employed in 
the theory of chemical reactions [T3 . [T5| and molecular 
motors [10, M, [igI [TtI and represent a significant sim- 
plification of the system dynamics. In this paper we be- 
gin with the continuous diffusion equation and systemat- 
ically transform to a discrete master equation for the case 
of overdamped Brownian motion on a multidimensional 
tilted periodic potential. 

There have been previous attempts to discretize the 
continuous diffusion equation for Brownian motion on a 
tilted periodic potential. In these treatments each cell of 
the periodic potential is identified with a discrete state 
and a master equation describes intercell hopping. It 
has been suggested that these discrete equations can be 
made consistent with the continuous theory by match- 
ing the dynamic structure factor [l^, [l^ , the eigenvalues 
20 | ] , or the ratio between backward and forward hopping 
211 . Alternatively, a discrete model can be derived by 
connecting each spatial region via appropriate absorbing 
boundary conditions and source terms ^22|], or by inte- 
grating the continuous probability density over a specific 
spatial region in the vicinity of each cell [l^, [H, [IJ] . In 
contrast to these treatments, we implement the classical 
analog of the tight-binding model of a quantum parti- 
cle in a periodic potential [25l - [28| . This is achieved by 



expanding the probability density of the continuous dif- 
fusion equation in a complete basis of Wannier states. In 
the limit of deep potential wells these states are localized 
in the minima of the periodic potential and provide a 
natural basis for the system. Using this basis we derive 
a discrete master equation that for deep potential wells 
is consistent with the physically-intuitive interpretation 
of hopping between wells. 

In the multidimensional case, Brownian motion on a 
tilted periodic potential is not tractable analytically in 
general and instead has been investigated numerically 
[29-32]. For non-separable potentials the tilt in one de- 
gree of freedom can induce a current in another and this 
coupling has important consequences in particle sorting 
experiments and molecular motors [8l.[l3ll30j. The tight- 
binding approach presented in this paper provides an an- 
alytic tool for treating multidimensional non-separable 
tilted periodic potentials. 

This paper is organized as follows. In section [TT] we 
introduce the Smoluchowski equation in our notation. 
In section IIIII we expand the probability density for the 
system in the Bloch eigenfunction basis for the periodic 
potential and consider the time evolution of the Bloch 
bands. In section ITVl we transform to the localized Wan- 
nier states of the periodic potential and derive a discrete 
master equation for the system. In section|V]we consider 
the limit of deep potential wells where the tilt dependence 
of the hopping rates can be determined analytically and 
hopping transitions occur predominantly between nearest 
neighbor wells. In section IVII we derive general expres- 
sions for the drift and diffusion. We conclude in section 

ivni 



II. SMOLUCHOWSKI EQUATION 

We consider Brownian motion of a particle in the over- 
damped limit of negligible inertia as described by the 
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Smoluchowski equation P, [T5| 
dP{r,t) 



dt 



(1) 



where P{r,t) is the probability density of finding the 
particle at position r at time t. The evolution operator 
for the Smoluchowski equation is 



V— — 

^ 7j dr 



e 



_9_ ay(r) 

9r, 9r, 



(2) 



where 1^(7") is the external potential, Q — ksT, ks is 
the Boltzmann constant, T is the temperature, and j is 
the coordinate index. The friction coefficient 7 may be 
different for each degree of freedom. 

The external potential V{r) is the tilted periodic po- 
tential 

V{r)^Vo{r)- f -r, (3) 

where Vo{r) is periodic with period a, i.e., 

Vo (r) ^Vo{r + a^r, )^Vo{r + a), (4) 

and the linear potential represents a constant macro- 
scopic force that drives the system out of thermal equi- 
librium and induces transport through the periodic po- 
tential. We impose periodic boundary conditions with 
a period given by Na where AT is a diagonal matrix of 
integers. We consider an infinite spatial extent by tak- 
ing the limit Njj — >■ 00, although a similar approach is 
possible for Njj finite. 



III. BLOCH EIGENFUNCTION EXPANSION 

The Smoluchowski equation ([1} has separable solutions 

Pk{r,t)^Mr)e-^''\ (5) 

where the eigenfunctions 4>kif) and eigenvalues satisfy 
the eigenequation 



C(j)k{r) = -XkMr). 



(6) 



Due to the periodicity of C, the eigenfunctions (j)k{r) can 
be chosen with the Bloch form [1, 28] 



(7) 



where Ua^k{i') has the same periodicity as the potential 
Vo{r), a is the band index, and k is continuou^ and 



restricted to the first Brillouin zone, i.e., 



TT , IT 
— <ki < — 



(8) 



^ For a finite spatial extent, Njj is finite and the label k is quan- 
tized. 



The evolution operator C can be transformed to the 
self-adjoint operator 1] 



n 
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where 



dV{r) 



2 9r? 



(9) 
(10) 

(11) 



The operator H is not Hermitian in general. However, 
denoting quantities for the f — case by a zero subscript 
(or superscript), the operator Tio is Hermitian and the 
untilted Bloch eigenfunctions </'°fc(^) form a complete 
orthonormal basis (see Appendix The untilted Bloch 
eigenfunctions satisfy the orthonormality relation 

dr e^°W/«</>°:,(r)e.fe'(r-) = <5„„.<5(fc - k'), (12) 

and the completeness relation 

eV-oW/e J2 f dk cj^l%{r)cl^lkir') - '5(r - r'), (13) 

a 

where the integral in Eq. ((T5)) is denoted by B to indicate 
that it is over a single Brillouin zone [see Eq. ([8])]. The 
eigenvalues A° ^ are real and exist in bands separated by 
forbidden band gaps (see Appendix [5| . 

Smoluchowski equation solutions satisfying periodic 
boundary conditions at infinity can be expanded with- 
out loss of generality in the untilted Bloch eigenfunction 
basis. The expanded probability density is 



where 



P{r.t)=Y, [ dk c„,fc(O0°,fc(r), 



(14) 



(15) 



Transforming the Smoluchowski equation ([T]), the coeffi- 
cients CaM{t) evolve according to 



dca^k{t) 
dt 

where 



/ ^ / dk fia,a' ,k,k' Ca\k' {t) : 



(16) 



dr e^°W/«0°\(r)/:e.'(r). (17) 



In the / = case, 

f^a,a',k,k' = ^^a,k^aa'S{k — fc'), 

and the evolution equation (jl6p becomes 
dclkit) 



dt 



^a.k^a,k 



it). 



(18) 



(19) 
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Equation ([19]) can be integrated directly to give 



(20) 



Equation ([SO]) shows that each eigenfunction coefficient 
c° ^ {t) decays in time with a rate given by its eigenvalue 
^ak — ^ (^^^ Appendix [BJ . The band gap between the 
lowest and first Bloch bands enables a separation of time 
scales so that for t ^ l/min(A°fc) all bands with q > 
are damped out and only the lowest a — band remains 
occupied. 

In the tilted periodic potential, 

fJ'a.a',k,k' = ->^aM^aa'6ik - k') - Va,a\k5{k - fe'), (21) 

and the evolution equation (jl6p becomes 

dCa,k{t) 



''^Q,feCQ,fc(t) —''^^Va,a',kCa',k{t)- (22) 



dt 



The first term on the right-hand side of Eq. is due to 
the untilted periodic potential Vo(r) and leads to decay 
of the coefficients Ca^k{i) at the rate A° ^ [see Eq. (^0]) ]. 
The second term is due to the linear potential and couples 
different Bloch bands of the untilted periodic potential. 
If the tilt / is sufficiently small that the band gaps in 
the decay rates A° are large compared to the coupling 

^a,a'^kj I.e., 



k\ 1 



(23) 



interband coupling can be neglected. In this weak-tilting 
regime, the evolution equation (|22p can be approximated 

by 



dt 



— ^[^a,k + '^a,a,k]Ca,k{t)- 



(24) 



Equation ([24]) can be integrated analytically and reduces 
to Eq. for / where Va^a,k — 0. 



IV. WANNIER BASIS 

In the case where the periodic potential Vb(r) creates 
deep potential wells that localize the particle, the untilted 
Bloch eigenfunctions (f>° ,^{r) are not a natural basis be- 
cause they are delocalized over the entire spatial extent 
of the system. It is more convenient to transform to the 
Wannier states that are localized around the potential 
wells of the periodic potential. The Wannier states for 
the eigenfunctions 0° ^(r) are 



Wa.n{r) ^ D dk fc(r)e 



-ik-An 



(25) 



where A is a diagonal matrix with Ajj = aj , the constant 
D is 



D 



(26) 



and for brevity we have omitted the zero superscript 
on the untilted Wannier states ([25]) . Taking the com- 
plex conjugate of the eigenequation ([6]) with / = 0, 
the untilted Bloch eigenfunctions are found to satisfy 
'/'Q*fc('') — 't^a -ki''') the Wannier states (|25)) are real. 
The untilted Wannier states form a complete orthonor- 
mal basis. They satisfy the orthonormality relation 

^Jdr e^°('-'/®w;„,„(r)u;„,,„,(r) = (27) 

and the completeness relation 

^e^oW/e ^ yjg.^r.ir)wgAr') - S^r - r'). (28) 

a, 71 

The probability density P{r, t) can be expanded in the 
untilted Wannier states as 



(29) 



■cMt) = / dr e^°W/^z«„,„(r)P(r,i). (30) 



where 



The coefficients Pa,n{t) are real and their amplitude can 
be interpreted as the probability that the Brownian par- 
ticle occupies the Wannier state Wa,n{'''). Transforming 
the Smoluchowski equation ([1]), the evolution of the sys- 
tem can be described by the master equation (l5l | 



dt 



where 



' a, a' ,n,n' 



I- [dr e^«W/^u;„,„(r)/:w;„,,„,(r).(32) 



D 



The coupling matrix defined by Eq. p2p is real and, 
recognising that the Wannier states satisfy 

Wa,n{r) ^ Wa,o{r - An), (33) 

we find that 

^a,cx',n,n' — f^a,a' ,n— n' ,0 ^a,a' ,0,n' — n- (^4) 

In the / = case, the coupling matrix becomes 
^°,a',r.,r.' = / (r )£oz«a' (r )(35) 

^a,n—n'^aa' 1 (36) 

where k° ^ are the Fourier components of the Cq eigen- 
values A° f^, i.e., 



K,°„ = -D I dk A° fcC^'^-^'- 



(37) 
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The master equation (|31]) becomes 



dt 



(38) 



and each band a evolves independently. The eigenvalues 
A° are real (see Appendix |A| and symmetric in k (due 
to the Bloch form), so by Eq. ([37)) the hopping rates k° ^ 
are real and symmetric in n, i.e., 

Equation (p9)) shows that the forward and backward rates 
of thermally activated hopping transitions between po- 
tential wells are equal. This is the condition of detailed 
balance [isj . 

In the tilted periodic potential, the coupling matrix is 



and the master equation ([31]) becomes 



dpa,n{t) 
dt 



(40) 



a' ,n' 



(42) 

The contribution A.a^a',n from the linear potential both 
modifies the intraband hopping rates and gives rise to 
coupling between the untilted Bloch bands. For strongly 
tilted periodic potentials the potential wells created by 
the periodic potential Vb(r) are substantially modified by 
the tilt and coupling between the untilted Bloch bands 
may be significant. The strong tiltiiig regime has been 
considered previously in both one [33l435j and two dimen- 
sions l30l. 



(see Appendix [C|) . i.e.. 



/ dk A„,fce''=-^'' 



The intraband hopping rates are real, i.e., 



(46) 



(47) 



but the linear potential breaks the symmetry so that ^ 
is imaginary in general (see Appendix \^ and Ka,n 
Ha~n- In the / = case, Eqs. ([43]) and ([46]) become 
Eqs. ([38]) and ([37]) . respectively. 

In the weak-tilting regime, coupling between the un- 
tilted Bloch bands is negligible and each band evolves 
independently. For t < 1/ A^ q , the system evolution is 
dominated by damping of the a ^ Bloch bands, as 
described in section Hill After this time the higher Bloch 
bands are barely occupied and the dynamics of the Brow- 
nian particle is dominated by intraband hopping within 
the a = band. This a = band preserves normaliza- 
tion under the evolution operator C (see Appendix [Pj) . 
Retaining only the stable a = band and dropping the 
a subscript, the evolution of the Brownian particle can 
be described by the single-band master equation 



dpnjt) 
dt 



Inverting Eq. (j46|) gives the eigenvalues 



n 

Taking A; = gives 



-ih-A.n 



-Aq.o, 



(48) 



(49) 



(50) 



A. Weak-Tilting Regime 

The Wannier basis is particularly useful when the force 
/ is not strong enough to significiantly modify the poten- 
tial wells of the periodic potential Vb(r) and interband 
coupling is negligible. This weak-tilting regime is defined 
by Eq. ([23]) and, while interband coupling is negligible, 
the linear potential modifies the intraband hopping rates. 
In the weak-tilting regime, the master equation <\A2\ can 
be well approximated by 



dpa,n{t) 
dt 



.'Pa,n' {t), 



where the intraband hopping rates are 

Ka,n = ^ J dr e^°^'^^/^Wa^ri{r)Cwa,o{r) 



(43) 

(44) 
(45) 



Neglecting interband coupling, the hopping rates Ka^n 
are the Fourier components of the C eigenvalues Xa.k 



and, using that Ao.o = (see Appendix ID]) . the master 
equation (|48p preserves normalization and can be cast 
into the more familiar form P, [isj 



dpnjt) 
dt 



— ^ ^ [^n—n'Pn'it) l^n' —nPn{ty\ 



(51) 



Equation ()51|) describes hopping between Wannier states 
within the a = Bloch band. With the hopping rates 
([46| . Eq. ([5T|) is a single-band discrete master equation 
describing Brownian motion on a weakly-tilted periodic 
potential and is the key result of this paper. 

The master equation (j5ip represents a significant sim- 
plification of the system dynamics and can be solved by 
transforming to the diagonal form 



dckjt) 
dt 



= -XkCk{t), 



where the eigenstates 



ik-Ar 



(52) 



(53) 



5 



are the expansion coefficients of the probability density in 
the untilted Bloch eigenfunction basis [see Eq. ([Tl)) ] and 
Afe are the eigenvalues of the evolution operator £ [20l |. 
in the weak-tilting regime where interband coupling is 
negligible. Equation ([52]) can be integrated analytically. 



B. Current 



The Smoluchowski equation ([T]) can be written as 
dP{r,t) 



dt 



-V-J{r,t), 



(54) 



where the coordinates of the current density are 



JAr,t)^ 

l3 



e 



d , dVf{r) 
drj drj 



Pir,t). (55) 



The current density can be expanded in the Wannier ba- 
sis as 



a,n 

with the expansion coefficients 

Ja,„(i)- J dr e^°W/%„,„(r)J(r,i). 
The current for the Brownian particle is 

1 = J dr J{r,t) = Y,jc..r,ma. 



(56) 



(57) 



(58) 



Equation ([58| is equal to the particle drift (see Section 
IVip when the current density vanishes at the system 
boundary [S]. 

In the weak-tilting regime retaining only the a = 
band, the current density expansion coefficients can be 
determined to be 

= -^Ep^'W;^ / '^^ ^r.(r)Ve^''M/Q«;„,(r) 



--Pn{t), 
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(59) 



where the vector division is taken elementwise. Taking 
the sum over n and using the Bloch form for the eigen- 
states with the ground state 4>QQ{r) oc e"^"'^''-'/^, the 
contribution to the current (j58p from the first term in 
Eq. ([5^ can be shown to vanish. Therefore, 



1 



(60) 



i.e., for weak-tilting and long times, the steady-state cur- 
rent is proportional to the tilt. 



V. TIGHT-BINDING LIMIT 

We consider the tight-binding limit where the Wannier 
states are strongly localized around a single period of 
the potential Vo{r) and the linear potential is sufficiently 
slowly varying that it changes negligibly across the spa- 
tial extent of any given Wannier state. With these con- 
ditions satisfied, it is possible to analytically derive the 
functional dependence of the hopping rates K^.n and the 
eigenvalues Xa,k on the force / and to further explore the 
dynamic behavior of the system without specific knowl- 
edge of the periodic potential. 

The validity criteria for the tight-binding limit can be 
determined formally as follows. The Wannier states are 
well localized if the periodic potential has deep potential 
wells, i.e., when the amplitude of the periodic potential is 
large compared to the thermal energy Q. With this con- 
dition satisfied the shape of the Wannier states (partic- 
ularly in the lowest Bloch band) depend predominantly 
on the curvature of the periodic potential at its extrema. 
Approximating the periodic potential in the vicinity of 
its extrema by a harmonic potential (see Appendix |E]) , 
the linear potential is a small perturbation near these 
extrema when 



I/, I «aj 



dr] 



(61) 



where Text the positions of extrema in the periodic po- 
tential. Equation (IFTI) ensures that the tilt varies slowly 
across the extent of the Wannier states and that coupling 
between Bloch bands is negligible (see Section IIV Ap . 

In general the coupling matrix aa.a' ,n,n' of Eq. (j32p 
can be expressed in terms of the self-adjoint operator H 
[see Eq. ^] as 



^ Q.a' .n.n' 



1 

D 



dr e^°«/«e-^W/2«u;„,„(r) 
xHe^W/2®«;„,,„,(r). (62) 



To evaluate Eq. (|62p in the tight-binding limit we insert 
the completeness relation (|28p and take the first term in 
the Taylor expansion of ±/ • r/2Q around the position 
An of the appropriate Wannier state, i.e.. 



1 

D 

This gives 

^ a. a' .n.n' 



dr e^°W/Qe±^-'-/2e^^ „(r.)^„, „,(r.) 



(63) 



-/•A(n'-n)/2e_ 

D 



xHe^°M/2e?«„,,„,(r). (64) 



The dominant contribution to the integral comes 
from the spatial regions around the extrema of the pe- 
riodic potential where the Wannier states are finite and 
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the potential is approximately harmonic |36j . In these 
regions the operator Ti can be approximated by Tio and 

^a.TL—n'^aa'^ (66) 



with 



fAn/2e 



(67) 



Equation (|66p shows that, in the tight-binding limit, the 
coupling matrix cra^a',n,n' is diagonal in the band index 
and the bands decouple in the same way that they decou- 
ple in the weak-tilting regime (see Section HV Ap . Equa- 
tion (|67p provides a simple functional dependence of the 
hopping rates on the tilt. This can be used to derive the 
ratio of forward to backward hopping rates, i.e.. 



(68) 



Equation ([68t is a generalized detailed balance condition 
for tilted periodic potentials d, [2l|, [2^ [23| ■ It is indepen- 
dent of the exact form of the periodic potential and the 
band index a. However, the Wannier states become in- 
creasingly delocalized with increasing band index so Eq. 
([68|) can be expected to become less accurate with in- 
creasing a. 

In the tight-binding limit, the tilt dependence of the 
Bloch eigenvalues Xa,k can be determined by inserting 
the tilt dependence (|57|) of the hopping rates into Eq. 
This yields 



A. 



^fAn/2e^-ikAn 



(69) 



In the case_of tight-binding, nearest-neighbor hopping 
dominates 27]. This means that higher-order hopping 
across multiple wells can be neglected and the summation 
in Eq. (|69p need only be extended over nearest neighbors. 

For long times t 3> l/A°oi the system dynamics is 
dominated by hopping within the a = band and can 
be described by the single-band master equation (|511) (see 
section HV Al) . In the tight-binding limit, A°q can be ap- 
proximated by the splitting of the two lowest eigenvalues 
of the harmonic approximation to the minima of the pe- 
riodic potential, i.e., A°o Z]j (^^^o(T')/9r||r=r„i„/7j, 
where rmin are the positions of the minima in the periodic 
potential. 

In the tight-binding limit, nearest neighbor hopping 
dominates and the summation in the master equation 
([5T|) need only be extended over nearest neighbors. The 
physical interpretation of the master equation can then 
be made more explicit by assuming that the Wannier 
states in the lowest band are the Gaussian harmonic oscil- 
lator states of the potential minima. With that assump- 
tion, the coefficients Pn(t) are non- negative according to 
Eq. (|5D|) . provided that the Wannier states are chosen 



positive. The coefficients Pn{t) can then be interpreted 
as the probability of the particle being localized in the 
nth potential well and the master equation ([ST|) describes 
particle hopping between nearest-neighbor wells. 



A. More General Potentials 

In the tight-binding limit, the physical arguments used 
to derive the tilt dependence of the hopping rates Ka^n 
can be used to better understand Brownian particle dy- 
namics in more general potentials. For example, consider 
the external potential 



V{r)^Voir) + Wir), 



(70) 



where Vo(t') is the periodic potential defined previously 
and the potential W{r) is slowly varying on the scale of 
a. In this case, the coupling matrix (Ja^oi',n,n' for the 
potential (|70)) can be approximated by taking the first 
term in the Taylor expansion of W(r) around the position 
An of the appropriate Wannier state, i.e., 

1 J dr e^«W/«e^W/2Q«;„,„(r)«;„,„, (r) 

^^w(A„)/2e^_,^_„ (71) 

Approximating the operator H by Hq yields 

W{An')/2e {An) /2e O r ^72^ 



and considering only nearest-neighbor hopping, Eq. ([72 
can be approximated as 



-VW{An)A{n-n')/2e 



I'Saa'- (73) 



To lowest order, VW{An) is independent of n and the 
dynamics are described using the tilted periodic potential 
already discussed. However, to higher orders Eq. ([75)1 
provides access to the dynamics for systems with external 
potentials of the form of Eq. (ffO)) . 



VI. MACROSCOPIC OBSERVABLES 

Macroscopic properties of the system can be calculated 
directly from the probability density P{r,t) or by using 
either the untilted Bloch eigenfunction basis or the Wan- 
nier state basis. We focus on the weak-tilting regime 
defined by Eq. ^ for long times {t > l/A?_o)- In this 
case the Wannier state basis is useful and the system is 
well described by the single-band master equation (j51l) . 
In the tight-binding limit the hopping rates are given by 

Eq. dnn). 
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Drift 



Diffusion 



The mean position can be calculated in the Wannier 
basis as 

{t) + 

(74) 



where 



Ra 



D 



dr rwa,oir) 



(75) 



In the weak-tilting regime retaining only the a = 
band, the time derivative of Eq. (|74p can be found by 
inserting the master equation (|51|) . This yields 

d{r){t) 



dt 



M ^ AriKr. 



(76) 



= -iAA VfcAfcl^^o (77) 

= A^VfcIm(Afc)|;,^o, (78) 
where J\f is the normalization. Integrating Eq. (|78p gives 

{r){t) = {rm+m VfcIm(Afe)|^^o , (79) 
and the drift is 



{rms 



lim ■ 



t 



A/" VfcIm(Afc)|fc^o . 



(80) 



Equation ((80|) shows that the drift is proportional to the 
gradient of the imaginary part of the lowest Bloch band 
eigenvalues [l^l- For the untilted case, VkA^ vanishes 
at fe = [see Eqs. ^ and ([39])]. Therefore, {ro){t) = 
(ro)(0) and the drift vanishes. This is a consequence of 
detailed balance for a periodic potential. 

In the tight-binding limit, the tilt dependence of the 
hopping rates and the eigenvalues Afc are known analyt- 
ically (see section \V\ and we find that 

d{r){t) 



dt 



1 



where the sum is taken over nearest neighbors. Consider- 
ing only one dimension labeled x, the net rate of nearest 
neighbour hopping is 

fx^x 



Kl — K_l 



-2ki sinh 



28 



(82) 



and the eigenvalues of C in the a = band are approxi- 
mately given by 

■f' Gj \ . ( Cli ^ , j'Ou ^ 



Afc^ = 4:Ki sin 



2 / V 2 
Differentiating Eq. (l83l) . the drift is 

(i)(i)st = 2Mn\a.x sinh 



+ 1- 



29 



fx^x 



(83) 



(84) 



Equation ([M]) is consistent with previous one- 
dimensional treatments [H, [33| . In two or more 
dimensions, Eq. (jSip shows that a tilt in one dimen- 
sion can induce current in another, provided that the 
potential is not separable [Tsj . 



The second moment can be calulated to be 
{r^){t) = J dr r^P{r,t) (85) 



where 



Sa = j dr r'^Wa,o{r). 



(87) 



In the weak-tilting regime retaining only the a = 
band, the time derivative of Eq. ((86|) can be found to be 



'^''^-(*)-i^+AA^(An)^.. 



dt 



N dt 

1 d{r)\t) 
N dt 
1 d{r)\t) 



NVlRe{\u)\ 



N dt 

It is convenient to define the variance 

var(t) = (r2)(i)-l(r)2(t), 

which evolves according to 
(ivar(i) 



fc=0 



dt 



= AAV^Re(A;,)|^^„. 
Integrating Eq. (|92|) gives 

var(t) = var(O) AAt ^l^<>^k)\k=o . 
and the diffusion coefficient 



(89) 
(90) 

(91) 

(92) 
(93) 
(94) 



Equation shows that the diffusion is proportional to 
the curvature of the real part of the lowest Bloch band 
eigenvalues (2^. In the tight-binding limit, the diffusion 
in one dimension is 



V = Nn\al cosh 



fx^x 

~2Q' 



(95) 



VII. CONCLUSION 

We have presented a systematic transformation from 
the continuous Smoluchowski equation to a discrete mas- 
ter equation using a Wannier basis expansion analogous 
to the tight-binding model of a quantum particle in a 
periodic potential. The discrete master equation can be 
interpreted in terms of intraband and interband hopping 
and, in the regime of weak tilting and long times, re- 
duces to a single-band master equation consistent with 
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the idea that the system relaxes rapidly within poten- 
tial wells and is then dominated by comparatively slow 
hopping transitions between wells [13] ■ The regime of 
validity is eq uivalent to Kramers' regime for a bistable 
potential j36l. Issj . 

In the tight-binding limit we have derived simple an- 
alytic expressions for the tilt dependence of the hopping 
rates and eigenvalues of the lowest Bloch band. These 
expressions have been used to determine the tilt depen- 
dence of the drift and diffusion in the long-time limit. 
When the periodic potential has a single dominant max- 
imum per period the hoppiiig rates can be approximated 
by Karmers' relation [ijj llSl, although our rates for 
f ^ are restricted to weak tilting so do not depend 
on the relative position of the potential maxima between 
minima ■ The tilt dependence of the ratio between 

forward and backward hopping rates takes the form of a 
generalized detailed balance condition that is consistent 
with nonequilibrium fluctuation theorems [40, 41]. 

In this paper we have demonstrated that the tight- 
binding approach is a valuable tool for deriving analytic 
solutions to Fokker-Planck equations with periodic po- 
tentials. This is particularly valuable for non-separable 
multidimensional potentials. The implications of this 
work for energy transfer in molecular motors have been 
highlighted elsewhere (T3| . 



Appendix A: Orthonormality of Eigenfunctions 

The self-adjoint operator H [see Eq. ([TU])] has the eigen- 
value problem 



where the eigenfunctions are 



.Vir)/2e 



Mr), 



(Al) 



(A2) 



and Afc and 4>k{f) are, respectively, the eigenvalues and 
eigenfunctions of the operator C [see Eq. Equa- 
tion (jA2|) shows that, when the eigenfunctions (pkif ) are 
chosen with periodic boundary conditions, the eigenfunc- 
tions ipkif) are not periodic. This means that for the 
tilted periodic potential the operator Ti is in general not 
Hermitian, the eigenvalues Xk may be imaginary, and the 
eigenfunctions ipkif) and (/)fc(r) do not necessarily form 
a complete orthonormal basis [1[ . 

For the untilted periodic potential with f — 0, the op- 
erator Co can be transformed to the self-adjoint operator 
Ho with eigenfunctions 



,Vo{r)/20 



(A3) 



Equation (|A3[) shows that both the eigenfunctions 4'ki''') 
and ipki''') simulaneously have periodic boundary 
conditions. For periodic boundary conditions with a pe- 
riod given by Na, the operator Tio is Hermitian. This 
means that the eigenvalues are real and exist in contin- 
uous bands separated by finite forbidden band gaps [28| . 



The eigenfunctions '0fc('") form a complete orthonormal 
basis . For an infinite spatial extent, the eigenfunctions 
-02 (r) satisfy the orthonormality relation 

dri:°*{r)^/j°,{r) = 6{k-k'), (A4) 

and the completeness relation 

dk^P°*ir)^lj°^ir')^S{r-r'). (A5) 



In the Co operator basis, the othonormality relation be- 
comes 

J dr e^°W/«0°*(r)0°,(r) = <5(fc - k'), (A6) 
and the completeness relation becomes 

e^^o(r)/e f 0O*(r)0O(r') = 5{r - r'). (A7) 



Appendix B: Positivity of Eigenvalues 

Following the methodology presented in Ref. the 
operator Co can be written as 



Co=Y. 



Q d 



.-Vo{r)/0 



d 



.Vo(r)/e 



(Bl) 



Using Eq. © with / = and Eq. (|X6)) . 

j dr e''°(^y'Ul*[r)Co4>l,{r) = -A°<5(fc - k') (B2) 



(B3) 



where the left-hand side of Eq. (|B2[) has been integrated 
by parts using Eq. (jBll) . Equation (jB3|) is less than or 
equal to zero so the eigenvalues A° are greater than or 
equal to zero. In particular, the ground state is 



with the eigenvalue Aq = 0. 



{B4) 



Appendix C: Eigenvalues for the Weakly-Tilted 
Periodic Potential 

The eigenfunctions 0Q,fe(r) of the tilted periodic po- 
tential can be expanded in the Wannier basis as 

(t>a,k{r)= E "^a^a'. ,k.n'Wa'.n'{r). (CI) 



Substituting expansion fjCip into Eq. ([6|, multiplying by 
gVb(r)/e^^^^^y,^^ and integrating gives 

-^Qj/c^Qja.fc.n — ^ ^ ^ a,a' ,k,n' ^a,a' ,n,n' • (^^) 
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Using the Bloch form for 0c<,fc(r), it can be shown that 



^a,a,k,n — ^ a,a' ,k,oe^'' ■ (C3) 



Substituting Eq. (jCSp into Eq. (|C2p and neglecting in- 
terband couphng yields 



Aa.fc = - ^ KQ,rie (C4) 

n 

which can be inverted to derive Eq. (|46|) . 



second term is proportional to Ao.a,n = ^'o.q.o and 

= 0, (D7) 

where we have used that the untilted Bloch ground state 
(t)l^o{r) oc e-^°^'^^'^ [see Eq. (HH)] is periodic at the 
boundary. This shows that the untilted a = Bloch 
band is stable when evolved by the operator C, even for 
a finite tilt. 



Appendix D: Normalization of the lowest untilted 
Bloch band 



Appendix E: Validity of the Tight-Binding Limit 



The normalization of the probability density P{r,t) 
can be calculated in the Wannier state basis as 

M{t) ^ Jdr P{r,t) = J2Po..r.ma, (Dl) 



where 



= y u;„,o(r). (D2) 



In the f — case, the normalization can be determined 
by taking the sum over n of Eq. ([55)1 and integrating to 
give 



(D3) 



Equation (|D3|) shows that the band a decays in time 
with a rate given by A° q > (see Appendix IB|) . Only 
the a = band with AfJ q = is stable. This is consistent 
with the separation of time-scales argument presented in 
section Hn] [see Eq. ([20| ]. 

For the tilted periodic potential, taking the time 
derivative of Eq. (jDl[) and inserting the master equa- 
tion (|42l) , the time derivative of the normalization can be 
found to be 



dt 



E 



^a,n'Pa,n 



.'Vo^'.nit) 



(D4) 

The contribution from the lowest a = band is given by 
the a = contribution to the sum in Eq. (jD4[) . The first 
term is due to the periodic potential and vanishes because 



-a;; = [see Eq. ^ and AppendixE]. The 



In the tight-binding limit, the shape of the Wannier 
states depends on the curvature of the periodic potential 
at its extrema and in the vicinity of these extrema the 
linear potential is a small perturbation. To derive the 
regime of validity for the tight-binding limit, the periodic 
potential in the vicinity of its extrema at r = roxt can be 
approximated by the harmonic form 



(El) 

'•'='-oxt 



The components of the potential L/°(r) of the Hermitian 
operator "Ho [see Eq. (|lip] can then be approximated by 



Uf{r) 



46 



(fj -re^tj? 



(E2) 



With the addition of the linear potential, the components 
of the potential U (r) become 



C/,(r) 



4Q ^oxtj j 2 



drf 



(E3) 



where the tilt shifts the extrema from r = Text to r = Text 
with 



' extj' 



' GX.tj 



d'Vo{r)/dr]l_^ 



(E4) 



The tight-binding limit corresponds to the regime where 
the shift I Text j — foxtjl is significantly smaller than the 
period aj of the periodic potential, i.e.. 



I/, I «aj 



drj 



(E5) 
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